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CN ! Abstract 

, A dynamical, non-Euclidean spacetime geometry in general relativity theory implies the possibility 

' of gravitational radiation. Here we explore novel methods of detecting such radiation from astrophysical 

^—i ' sources by means of matter- wave interferometers (MIGOs), using atomic beams emanating from super- 

I . sonic atomic sources that are further cooled and collimated by means of optical molasses. While the 

y—^ ' sensitivities of such MIGOs compare favorably with LIGO and LISA, the sizes of MIGOs can be orders 

of magnitude smaller, and their bandwidths wider. Using a pedagogical approach, we place this problem 
■ into the broader context of problems at the intersection of quantum mechanics with general relativity. 

> ■ 

0\ . 1 Introduction 
O ■ 

. In this contribution to the proceedings of the 2003 Gargnano workshop on 'Mysteries, Puzzles, and Paradoxes 

' of Quantum Mechanics', let us begin with a brief review of those basic notions of differential geometry that 

, will be important for the understanding of the interaction between gravitational radiation and quantum 

' matter. Specifically, we shall examine here matter-wave interferometry as a means of detecting gravitational 

radiation from astrophysical sources. In a pedagogical manner, we shall place this problem into the broader 
^ ' context of some of the fascinating problems at the intersection of quantum mechanics with general relativity, 

{^j^), One of Euclid's axioms for plane geometry is the parallel-lines axiom, i.e., that parallel lines never meet, 

^ ' no matter how far these lines are extended in either direction. This axiom is illustrated by figure ^a). We 

• 1— I . shall presently extend this axiom of Euclidean geometry to a statement about worldlines of test particles in 

' flat spacetimes in general relativity. 

^ , Gauss discovered non-Euclidean geometry after a deliberate omission of this Euclidean axiom ^ |2] . 

- - ' He found one example of this geometry in the surface of a sphere, where there exists a constant positive 

curvature everywhere. This Gaussian geometry is illustrated in figure ^b). A violation of the Euclidean 
axiom can be illustrated by extending northward two initially parallel 'lines of longitude' starting at points O 
and P at the equator. These lines of longitude are geodesies, i.e., the 'straightest' possible lines between two 
points, which minimize the intervening distance between them. When these geodesies are extended north, 
they will converge onto, and eventually intersect at, the north pole N, thus violating Euclid's parallel-line 
axiom. 

Riemann generalized Gauss's notion of non-Euclidean geometry to that of differential geometry for 
arbitrarily curved manifolds 0. As a measure of local curvature, he introduced the operational method of 
the parallel transport of a vector around an infinitesimal circuit, i.e., a tiny closed curve on a differentiable 
manifold. After the parallel transport is completed, there results a rotation of the final direction of the 
parallel-transported vector by an angle with respect to its initial direction, as illustrated in part (b) of 
figure ni This rotation angle </) is thus a measure of the local curvature of the manifold. 

To quantify this angle Riemann introduced his curvature tensor R\ji^, which we define following 
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Figure 1: (a) Euclid's parallel-lines axiom in plane geometry states that two straight, parallel lines, such 
as the ones extended through the nearby points O and P, never meet at infinity. Gauss discovered a non- 
Euclidean geometry which is illustrated in (b). Two parallel geodesies on a sphere, represented by two lines 
of longitude emanating from points O and P on the equator of the sphere, converge onto the north pole 
N. Parallel transport of a vector around the spherical triangle ONP, and back to O again, illustrates the 
Gauss-Bonnet theorem, which states that the sum of the three interior angles of a triangle is 180° plus the 
solid angle subtended by this triangle with respect to the centre of the sphere. Here (j) is the rotation angle 
of the parallel-transported vector after its completion of the circuit. In Berry's phase, an 'anholonomy' is 
the quantum phase factor exp{ims(t>) that a spin with component mg normal to the surface would pick up 
after adiabatic transport around the circuit OPNO in the field of a magnetic monopole located at the centre 
of the sphere. 



Landau and Lifschitz's approach as 

- -iRljkSf'^l (1) 

where the Latin indices i,j,k,l, which run from 1 to 2 to 3, represent the three spatial dimensions, ^/ are 
components of the initial vector, and S^i are the changes of the final vector components relative to the 
initial components, after the process of parallel transport of the vector around a closed circuit — with a 
infinitesimal area given by Sf^'^ = Sx^Sx'' — is completed (Einstein's summation convention is used here for 
repeated indices). 

Einstein generalized Riemann's notion of curved space to curved spacetime by the generalization of the 
curvature of space to that of the curvature of spacetime, so that 

where the Greek indices A, /x, i^, n now run from (representing time), to 1, 2, and 3 (representing space). He 
then applied the geometrical concepts of Gauss and Riemann to describe gravity: Matter acts on geometry 
by curving spacetime, and geometry acts on matter by determining the paths of test particles Flat 
spacetimes still obey Euclid's parallel-lines axiom, now generalized to the case of two straight and parallel 
worldlines, i.e., trajectories of noninteracting particles in space and time, for two small, nearby free objects 
at rest with respect to each other, e.g., in outer space, located, at time t = 0, in their proper frames at the 
two points O and P, as illustrated in figurej^Ja). These two straight and parallel worldlines are also implied 
by the rectilinear, uniform motion of free objects obeying Newton's first law of motion. 

In this space-time diagram, time, represent by the vertical axis, ascends vertically, and space is repre- 
sented in one dimension by the horizontal axis. Thus one of the objects, say the observer O, is represented 
by a vertically ascending, straight worldline, and a nearby object P is represented by the nearby, parallel 
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worldline. If spacetime were strictly flat (as shown in figure|2Ia)), these two straight and parallel worldlines 
would never meet. However, like the great circles on a sphere, or geodesies on a pseudosphere, in curved 
spacetimes these lines can either converge, or diverge from each other, respectively, depending on the sign 
of the components of the Riemann curvature tensor. An observer fixed on O would see the object P either 
come towards him, or go away from him, and would interpret this motion as being due to an effective 'tidal' 
force acting on P (see Chapter 1 of ^ for a description of this). In figure |2Ib) we have drawn the specific 
case where components of the local Riemann curvature tensor are negative, and the worldlines diverge from 
one another. Since both the observer and the object are simply following their own worldlines — which are 
geodesies in the absence of all other forces on both the observer and the object — the equivalence principle 
still holds, and consequently this force is proportional to the mass of the object. The resulting motion of 
P as observed by O is therefore independent of the mass, composition, or thermodynamic state of the test 
object at P, and herein lies the geometrical meaning of the equivalence principle. Moreover, when the local 
Riemann curvature tensor between O and P can be approximated as a constant, this force, like the action 
of the Moon's gravity on the tides, increases linearly with the distance between O and P (see figureEfb)), 
and indeed, the use of the name 'tidal' for this effective force comes from the tidal force on the oceans by 
the Moon's gravitational field. 

It is important to note that while the worldline illustrations in figure|21are a visualization of geometry, and 
its connection to gravity and gravitational tidal forces, they are only a convenience, and do not represent 
any frame that can be achieved physically. As drawn, the perspectives are that of an Observer who is 
removed from the spacetime, and standing outside of it looking in. This is unphysical. Every experimental 
measurement is made through an experimental apparatus — which may be as simple a device as a pair eyes 
used to see the motion of an object P a small distance away — that is, by necessity, a physical object. As 
such, the apparatus must lie along a worldline in the spacetime of the universe, and cannot be removed from 
it. No experimental measurements can be done by a fictitious observer 'suspended' outside of the universe. 

Of particular interest to this paper is the case when a 'gravitational wave' (GW for short) is present in 
the spacetime. In this case the worldlines shown in figure [2Ib) will periodically converge and diverge from 
each other, due to the fact that the curvature of spacetime is alternating between positive and negative 
values as the GW passes over the observer O and object P. Such distance changes between the two objects 
can be measured by sending a light beam out from O, and reflecting it by means of a mirror at P back to 
O. The distance between the two objects, measured in this way, will also change periodically in time when 
the GW passes over them. Once again the observer O would interpret the motion of the object P as due to 
an effective 'tidal' force which, in the limit of long-wavelength GWs, also grows linearly with the distance 
separating O and P (see Chapter 35 of 0] for a more detailed description). 

Because of the experiment by Pound and Rebka |S] to measure the gravitational redshift, along with the 
experiment by CoUela, Overhauser and Werner (COW) [HI to demonstrate through neutron interferometry 
that the gravitational potential of the Earth induces a quantum phase shift, the general relativistic effect 
most familiar to the Atomic, Molecular, and Optical (AMO) community is the gravitational redshift; GWs 
and their interaction with matter are not nearly as well known. There is thus an unfortunate tendency 
among the community to try to understand all general relativistic effects, including those due to GWs, in 
terms of the redshift. This cannot be done. If the redshift were used in an attempt to understand the physics 
of GWs, it would be like trying to understand electromagnetic (EM) waves using the scalar potential only; 
most of the essential physics would be lost. 

Because of the AMO community's unfamiliarity, we present here a review of GWs. This review is 
necessarily brief; the literature on GWs, their generation, propagation and interaction with matter, is vast 
(see ^ and [7|), and a detailed presentation on the topic cannot be done within this paper. We will 
focus here on the physics that underlie GWs, and will delay the presentation of a mathematical review to 
Appendix A. 

The existence of GWs is one of the most fundamental predictions of general relativity. Unlike the 
Newtonian gravitational potential, which can be obtained from general relativity by taking the appropriate 
non-relativistic and weak gravity limits, GWs do not have an analog in Newtonian gravity. This is because 
like electrostatics, Newtonian gravity is a static theory that cannot encompass wave propagation. Indeed, 
the analysis of the gravitational redshift, as measured in the Pound- Rebka experiment, could be done using 
only special relativity and the equivalence principle; it does not need the full structure of general relativity. 
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(a) (b) 

Figure 2: Euclid's parallel-line axiom applied to the flat spacetime shown in (a) implies that the worldlines 
of two small objects O and P at rest with respect to each other at t = 0, remain parallel to each other 
and never meet. In the curved spacetimcs shown in (b), however, the two worldlines can either converge or 
diverge, or, as in the case of GWs, can oscillate. The acceleration of an object seen by the observer O close 
to P is proportional to the local Riemann curvature tensor, and is also proportional to the separation 
between O and P, as is the case in all such tidal effects. 



Just as it took the discovery of Ampere's and Faraday's Laws to lead Maxwell to predict the existence of 
EM waves, it took Einstein's generalization of Newtonian gravity to predict the existence of GWs, and to 
describe their properties. 

Many, but not all, of the underlying properties of GWs in the linear approximation can be understood |H] 
in analogy with EM waves (see Appendix A). Like EM waves, GWs have two physical polarizations — the + 
polarization and the x polarization — but unlike an EM wave, which is a spin 1 field and can be represented 
by a vector potential a GW is instead a spin 2 field and must be represented by a second-rank tensor 
h^y. This tensor is analogous to a local strain induced on spacetime by a GW, and is a measure of the size 
of the 'ripples' in spacetime caused by the passage of the GW. Like A^, h^i, is 'gauge' field, with the choice 
gauge being a choice of coordinate systems, rather than a choice in a U{1) phase factor. The gauge invariant 
object corresponding to the field strength for EM is the Riemann curvature tensor R^^^, constructed 
from second derivatives of h^iy. While it is possible to choose a gauge where /iqo 7^ — the term that would 
generate gravitational-redshift-like effects — this would not be the minimal gauge for the GW. Making this 
gauge choice would be analogous in EM to making a gauge choice where / instead of the usual 
radiation gauge: ^0 = and V • ^ = 0. Moreover, since the potential /iqo can describe at most one of the 
two physical degrees of freedom of the GW, attempting to understand the effects of GWs on matter using 
this /loo through a gravitational redshift argument would miss much of the underlying physics. Both degrees 
of freedom for the GW must be taken into account, and this is most easily done by choosing a minimal 
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Figure 3: The force- field lines for the + polarization (a) and x polarization (b) of a GW (from jiZ2i). These 
patterns are the snapshots of anisotropies of space, or 'strains', induced by the GW in the spacetime manifold, 
taken at i = by an observer located at the centre, i.e., snapshots of the patterns of motion of an ensemble 
of freely- falling, noninteracting test particles distributed uniformly over the (x, y) plane. These time- varying 
patterns of anisotropics cause time-varying phase shifts to appear both for light waves in Michclson-type 
interferometers, such as in LIGO, and also for matter waves in Mach-Zchnder-type interferometers, such as 
in MIGO. 

gauge. This gauge — which is analogous to the radiation gauge for EM waves — is the transverse-traceless 
(TT) gauge: = 0, V^/i^i, = 0, and for GWs propagating on flat ('Minkowski') spacetime, /iq^ = 0. 
In particular, the gravitational potential vanishes since /iqo — 0. There can be no gravitational-redshift- 
like effects; only the spatial 'strain' components hij are non- vanishing. Like A^, in this gauge hij can be 

represented as plane waves with two polarization vectors e^^"^ and , shown graphically on figure |3| 

One aspect of GWs that cannot be understood in analogy with EM waves is the interaction of GWs 
with matter. The reason for this goes beyond the obvious differences in the tensorial ranks of h^^ and A^. 
Unlike EM waves, GWs, like all gravitational effects, cannot be screened. As mentioned above, all physical 
observers must be part of spacetime, and cannot be removed from it. When a GW passes through a physical 
system, it acts on all parts of the system, including the observer, and consequently, only differences in the 
motion between the observer (O in figurel^J and the observed (P in figure^) can be measured; the 'absolute' 
motion of either one cannot (see discussion in Chapter 35 of 111). 

If only relative motion can be measured when a GW is present, what coordinate system should be chosen? 
A description in words of this coordinate system has been given by Thorne 0, and has been outlined by 
Synge and de Felice and Clark 11 for general spacetimes. More recently, Speliotopoulos and Chiao 
|12| has explicitly constructed this coordinate system for general spacetimes in the limit of linearized gravity. 
Its construction is based on the following physical constraints I3j : Every physical particle travels along a 
worldline in the spacetime of the universe. Every measurement of the physical properties of the test particle 
by an observer must be done using an experimental apparatus. Before the observer can take measurements 
with this apparatus, he must first choose a local orthonormal coordinate system. In curved spacetimes, 
this involves the construction of a local orthonormal coordinate system (a tetrad frame) |10| . Naturally, 
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this coordinate system will be fixed, say, to the centre of mass of his experimental apparatus, and will thus 
propagate in time along the worldline of the apparatus as well. The observer uses the coordinate time of the 
physical apparatus to measure time, which, because he or she will not be moving relative to the apparatus, 
is also his or her proper time. Thus the time axis of the coordinate system he or she has chosen will always 
be tangential to his or her own worldline. The position — which can be of finite extent — of the test particle is 
measured with respect to an origin fixed on the apparatus, and is the shortest distance between this origin 
and the particle. However, because the apparatus travels along its worldline, the origin of the coordinate 
system will also travel along a worldline in the spacetime. When the rate of change of the position of the 
particle is measured at two successive times, the relative four-velocity of the particle with respect to the 
apparatus will naturally be obtained (see ^21 for this construction). 

As described, it would seem that the observer simply constructs his or her usual laboratory frame, and it 
is difficult to see the novelty of this construction. The construction is done, however, in a curved spacetime, 
and as a result, in this coordinate system the observer O sees the object P in figure |2Ib) undergoing an 
acceleration dependent upon the local Riemann curvature tensor of the spacetime. When no other forces 
are acting on the object, it is straightforward to show that when P is close to O, this acceleration is 

dt^ - ^0.0, a; - 2 ^ ' y-^) 

where i?ojOj are components of the Riemann curvature tensor. The second equality holds for the special 
case when long-wavelength GWs (see Appendix A) — which causes the deviations hij of the metric from 
flat spacetime — are present. This equation is obtained by taking the usual geodesic equation for the object 
P, subtracting from it the geodesic equation for the observer O, and expanding the result to first order in 
the distance separating O and P shown in figure [^Jb). Equation ijSJl is the geodesic deviation equation of 
motion, and describes the difference in the motion between two nearby geodesies. Note that the acceleration 
in Eq. Q is proportional to Roioj, and is a gauge-invariant quantity. Consequently, although the explicit 
dependence of this acceleration on hij depends on the choice of the TT gauge for the GW, the dynamics (or 
acceleration) of the object described by Eq. Q will not, in the end, depend on this gauge (or coordinate) 
choice. 

Turning now to geometry and quantum systems. Berry pointed out that non-Euclidean geometry can 
often be encountered in quantum mechanics jT5j as well as classical mechanics. For example, a spin-half 
particle can be constrained by a strong magnetic field to always remain parallel to the direction of the field. 
When the direction of the magnet that produces this field is slowly rotated through space around a closed 
conical circuit, there results a rotation of the spin through the angle after the completion of a circuit. 
After this round trip, the spin picks up a quantum phase factor which is a consequence of the non-Euclidean 
geometry of the Bloch sphere (see figure |3b)). 

This phase factor can be explicitly seen through the Gauss-Bonnet theorem, which states that (p — fl, 
where is the solid angle subtended by the circuit with respect to the centre of sphere. There therefore 
results a quantum phase shift A$ experienced by the rotated spin after this round trip, where A<I> = 
niscj} = rrigQ for the case where the spin stays adiabatically in its eigenstate |ms). After the circuit is 
completed, the spin thus picks up a phase factor of exp {ims(j)) — exp [imris^) that can be measured in 
quantum intcrfcromctry. This is Berry's 'geometric phase factor', or 'anholonomy', which is similar to the 
Aharonov-Bohm phase factor picked up by an electron after it traverses a circuit on the surface of a sphere 
surrounding a Dirac monopole. 

Irrespective of the great successes of difi^erential geometric concepts in classical physics, since all ob- 
jects at a fundamental level obey quantum mechanics, the notion of 'worldlines', or of 'geodesies', or of 
'trajectories' in general, is problematic because of the uncertainty principle. As Bohr has taught us, we 
must fundamentally abandon the notion of classical trajectories in quantum mechanics. Nevertheless, from 
Berry's work it is well known that geometric concepts are fruitful in quantum mechanics as well. The crucial 
question then arises: How does one operationally measure spacetime curvature in a quantum world? As a 
corollary, how does one know whether a spacetime is precisely flat or not? The answer: By means of the 
interference of quantum test particles using matter- wave interferometry, as these test particles move through 
spacetime. These test particles are quantum in the sense that in the interferometer, we cannot know, even 
in principle, which path, in the quasi-classical limit, a given particle actually took. Hence we are strongly 
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Figure 4: SQUIDs as detectors of long-range fields. 'X's denote Josephson junctions, (a) A superconducting 
SQUID is a sensitive detector of magnetic fields through the Aharonov-Bohm effect, (b) A neutral SQUID 
(e.g., using superfluid helium) could also in principle detect gravitational radiation through a gravitational 
version of the Aharonov-Bohm effect. As a GW with + polarization passes over the interferometer propa- 
gating from left to right, the left and right sides of the SQUID, which are spaced by half a wavelength of 
a GW, causes tidal forces that reverse in sign from left to right. These forces induce a time- varying phase 
shift, which can be detected by means of the Josephson junction. 

motivated to examine the problem of matter-wave interference as a means of quantifying the problem of the 
quantum measurement of the curvature of time- varying spacetime geometries, and in particular, the specific 
case of those associated with gravitational radiation. This is closely related to the more general question: 
How does gravitational radiation interact with quantum matter? As we shall see below, this question 
already arises in an interesting way at the level of nonrelativistic quantum matter interacting with weak, 
linearized GWs. This problem is a linear one both in the matter and the field sectors of the theory, and 
should therefore be calculable using linear response theory. 

One of the authors (RYC) started thinking about this problem two decades ago when he asked 
the question: If the superconducting quantum interference device (SQUID), an early form of a matter- wave 
interferometer, can very sensitively measure the Aharonov-Bohm phase, what would happen if one were to 
replace the charged particles inside the SQUID by neutral particles? Is there some magnetic-like field leading 
to an analog of the Aharonov-Bohm effect, so that the neutral version of the SQUID can also sensitively 
measure this field? At that time, the answer seemed to be that the Lense-Thirring field of general relativity 
was this magnetic-like field. A natural extension of this question was: If one could sensitively measure 
stationary gravitational fields by means of quantum interference, could one also sensitively measure time- 
varying gravitational fields, i.e., gravitational radiation? One answer: The quadrupole (or tidal) symmetry 
of gravitational radiation would require that the neutral SQUID be twisted into a figure-8 shape before it 
could detect such radiation (see figure E^b)). 

Since the neutral SQUID, viewed as a macroscopic quantum system, would become, for weak gravita- 
tional fields, linear and reciprocal in its response to GWs, the question naturally arose: Could this figure-8 
antenna also generate, as well as detect, GWs? Recently, RYC suggested a simplification of these earlier 
ideas, which eliminated the hard-to-make Josephson junctions |17j . Quantum systems such as the quantum 
Hall fluid, a ferromagnetic material that responds linearly and reciprocally to external perturbations, could 
lead to a kind of 'quantum transducer action', in which conversion of GWs into EM waves, and vice versa, 
could in principle occur. This coupling originates from an extension of the usual minimal coupling rule to 
include the coupling of spin to curved spacetime |17j . 
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For example, the spin of the electron in a quantum Hall fluid can be viewed as a handle, by which it 
can be twisted around periodically in direction by a GW, thus causing it to pick up a periodic Berry phase. 
(The electron spin, which is a gyroscope-like object, locally undergoes parallel transport in the presence of 
a GW.) The Berry phase in turn induces macroscopic quantum flows in this charged quantum fluid which 
radiate EM waves. This macroscopic quantum process, in which a GW generates an EM wave, is both 
linear and reciprocal. Because of time-reversal symmetry, the time-reversed process in which an EM wave 
generates a GW, must also occur with equal power conversion efflciency; this efficiency can be at most 
unity . This reciprocity principle suggests the possibility of a Hertz-like experiment for both generating 
and detecting high-frequency GWs, in which, for example, microwave EM radiation is first converted into 
a GW by one sample of the quantum fluid, and then the generated microwave GW is back-converted into 
EM microwaves by a second sample of the same quantum fluid. (Faraday cages would prevent the usual 
EM coupling between the generation and detection parts of the apparatus.) A first attempt at this type 
of Hertz-like experiment was performed using YBCO, a superconductor instead of a quantum Hall fluid, at 
liquid nitrogen temperature as the quantum sample |2()j . An upper limit on the power conversion efficiency 
for the quantum transducer action of the YBCO sample was placed at 15 parts per million. YBCO is a zero- 
spin superconductor, however. Better choices might have been non-spin-zero materials like the spin-triplet 
superconductors, such as Sr2Ru04 [2], or the ferromagnetic superconductors, such as URhGe j22j. 

2 Matter-wave interferometry 

Our primary motivation for studying GWs and their interaction with matter on the quantum level is 
based on the observation that matter-wave interferometers can be very sensitive detectors of GWs from 
astrophysical sources. The rest of this paper will be focus on demonstrating this sensitivity, and exploring 
the suitability of using matter- wave interferometers to construct MIGO, the Matter- wave /nterferometric 
Gravitational-wave Observatory. 

The gain in sensitivity expected in using matter-wave interferometry instead of laser interferometry in 
detecting GWs can be seen from the following argument: Roughly speaking, since an atom with mass m 
'weighs' much more than a photon with frequency w-y, all other things being equal, a matter- wave-based 
Sagnac interferometer will be mc^/hojj ~ 10^° times more sensitive as an inertial sensor of Earth's rotation 
than a light- wave-based interferometer [23 • This fact is well known , and the advantage of an atom-based 
Sagnac interferometer over a laser gyroscope has recently been shown experimentally [JSlEni- Although 
other important factors will also play a role, the inherent advantage of matter-wave interferometers over 
light-wave interferometers is expected to be pervasive, and not restricted only to the Sagnac interferometer. 

Matter-wave interferometry is based on the particle-wave duality of quantum mechanics. This duality 
states that any massive object, such as an atom, can, under certain circumstances, behave like a particle, but 
can, under different circumstances, behave like a wave with a deBroglie wavelength X^b = 2Trh/mv, where 
V is the speed of a nonrelativistic {v << c) object. Like the phase of light, the quantum phase of an atom 
can be used to construct an interferometer using atoms, and we shall show that very sensitive observatories 
for GWs, which are expected to be many times smaller than corresponding light- wave- interferometer-based 
GW observatories such as LIGO (Laser Interferometer Gravitational-wave Observatory) and LISA (Laser 
Interferometer Space Antenna), can be constructed with matter-based interferometers. 

We shall argue here that the rapid technological advances in AMO physics that have occurred since 
the pioneering .27 work on atom interferometry have now made it possible to construct MIGO. These 
technologies did not exist in the 1970s and 80s when LIGO was first conceived and developed (see [7j for 
a history of LIGO); the inherent advantage of matter- wave interferometers over light-wave ones — suggested 
by the factor mc^/huiy ^ 10^° — either was not widely known or was not appreciated. Yet due in great 
part to this factor, MIGO can not only be many times smaller than either LIGO or LISA, but can also 
expand the frequency response range of current GW observatories to presently inaccessible regions, and can 
substantially extend the observational reach of GW observatories as well. It will also be possible to measure 
with MIGO other general relativistic effects that are not accessible to LIGO or LISA. 

The study of matter-wave interferometry in connection with GWs was started by Linet and Tourrenc 
[2H!- A few years later, Stodolsky also studied matter- wave interferometry and its use in measuring general 
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relativistic effects, including the detection of GWs. Both of these studies were done on the quantum me- 
chanical level, and used as their underlying basis the geodesic equation of motion, not the geodesic deviation 
equation of motion, to calculate the phase shift of individual particles in both stationary spacetimes, and 
spacetimes with GWs. These works were complementary to those by Anandan |3UII31[I5^ . who focused more 
on matter-wave interferometry in stationary spacetimes and the Sagnac effect. He also used the quantum 
mechanical framework based on the geodesic equation. More recently, statements of the impracticality of 
using atom interferometry to detect GWs have been made by Kasevich and coworkers [21] . Studies of GW 
antennas using quantum condensed matter systems such as superfiuids have also been done, and were started 
by Anandan and Chiao \6'6\ . The extension of the quantum- mechanical-based analysis to a quantum- field- 
theoretic approach for spin-zero particles was done by Cai and Papini |34|. Later, Borde and coworkers, in 
a series of papers I36L 157] , revisited the question of matter- wave interferometry and the measurement of 
gravitational effects — both stationary and nonstationary — using also quantum-field-theoretic methods, but 
now with the Dirac equation on curved spacetimes for a spin- half particle as their starting point. All of these 
quantum-field-theoretic approaches nevertheless resulted in a phase shift that is similar to Linet-Tourrenc's 
and to Stodolsky's. 

What has not been well appreciated in the previous work in matter-wave interferometry and its use in 
detecting GWs is that gravity, unlike the other forces of nature, cannot be screened. Thus, when a GW passes 
through a system, say a matter- wave interferometer, the GW acts on all parts of the system. Consequently, 
only differences between particles can be measured |^ 112). and in the long- wavelength limit for GWs, the 
motion of particles is described by the geodesic deviation equation Eq. |(2Jl, and not by the geodesic equation 
on which previous analyses have been based. This fundamental feature of GWs is well known in the general 
relativity community, and is explicitly exploited in the design of both LIGO and LISA. In this paper, we 
will follow LIGO and the standard general-relativistic approach of using the geodesic deviation equation as 
a starting point in our analysis of matter- wave interferometry in connection with GWs. As a consequence, 
we arrive at an a phase shift for MIGO that is very much different than those calculated before based on the 
geodesic equation of motion. A detailed comparison of the approaches based on the geodesic equation, such 
as those listed in the above, versus our approach based on the geodesic deviation equation will be presented 
elsewhere |38| : we shall only present a general critique of the geodesic-equation approaches here. 

Unlike LIGO, which uses laser interferometry to very accurately measure slight shifts in the position of 
a classical test mass — the end mirrors of the interferometer — when a GW passes by, MIGO measures slight 
shifts in the velocity of atoms caused by the GW. These velocity shifts in turn produces small changes in 
the deBroglie wavelength of the atom — a quantum test mass — and thus in its quantum phase. Like laser- 
based interferometers, these phase shifts can then be detected using interferometry, and we shall show that, 
assuming equal shot-noise limits, 



where Lmigo is the effective size of MIGO, -B = 75 is the number round trips the light beam takes within 
the Fabry-Perot interferometer placed in each arm of LIGO, cu^ is the frequency of LIGO's lasers, m is the 
mass of the atom used in MIGO, / is the frequency of the GW, and Ll/go is the physical length of one of 
LIGO's arms. At 125 Hz, where LIGO I, the current LIGO configuration, is most sensitive, Lmigo ~ 1-5 
m as compared to Lligo = 4 km if a ^^'^Cs atom is used Similarly, a comparison of a space-based 

MIGO with LISA leads to a MIGO configuration that is at least ten thousand times smaller than LISA, if 
they both have the same shot-noise sensitivity. 

It should be emphasized that at a fundamental level MIGO measures the local Riemann curvature tensor 
710 matter what its source is. This can most clearly be seen by drawing in spacetime the paths of an atom 
passing through an interferometer as was done in j42| : the two possible paths of the atom naturally forms a 
closed path in spacetime, similar to the ones is used to define the Riemann curvature tensor (see figuremb)). 
It is thus not surprising that like the well-known Berry's phase JHl and the Aharonov-Bohm effect for GWs 
|12| . the Riemann curvature tensor encircled by the closed path induces a net quantum phase shift. However, 
for GWs we shall see that the additional action of the mirrors in MIGO on the atom produces a larger phase 
shift than that caused by the curvature tensor alone. Thus, in addition to detecting GWs, MIGO, unlike 
LIGO, will also be able to measure directly the local Riemann curvature tensor from stationary sources such 




(4) 
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as the Earth, Moon, and Sun for the first time, and it may also be possible to measure the Lense-Thirring 
field of the Earth as well. 

The main objective of this paper is to present the following: The concept of MIGO along with its un- 
derlying theoretical framework; an outline of the design of two different types of MIGO configurations; a 
calculation of the expected sensitivities of these configurations; an assessment of their potential as gravita- 
tional wave observatories. We shall also argue as to the feasibility of constructing a MIGO. To this end, 
using our analysis of the underlying physics of MIGO, we have estimated the specifications of an interfer- 
ometer that would be capable of measuring GWs, and have outlined the technologies that could be used 
in reaching them. What has not been done here is a complete analysis of the systematic errors of MIGO. 
Since such an analysis would depend crucially on the precise structure of the interferometer, it would not 
be fruitful to do this without at least a small-scale MIGO in hand. We have, however, estimated the funda- 
mental limitations to MIGO's sensitivity due to thermal fluctuations of the mirror. Moreover, as we shall 
see below, the various technologies needed to construct MIGO have already been separately demonstrated 
in various atom diffraction and interferometry experiments since the early pioneering work in the 1990s. 
What has not been done is the integration of these components into a complete atom interferometer that 
has a design and a sensitivity to detect, measure, and observe GWs. 

It is unfortunate that misconceptions of atom interferometry, on the one hand, and laser gravitational 
wave detection schemes, on the other hand, are prevalent among both the LIGO and the AMO communities. 
Among the AMO community, on the one hand, it is often believed that it is the light beams in the arms 
of LIGO's interferometer that is being acted upon by the GW; the mirrors plays little role. Precisely the 
opposite is true. In Thome's 'proper reference frame', it is the action of the GWs on the end mirrors of 
the interferometer, which are suspended vertically with piano wires, that is being measured; light is used 
simply as a means to measure the shifts in the mirror's position. In addition, the AMO community often 
tries to understand the properties of GWs and their interaction with matter through an analogy with the 
gravitational redshift. This also is too naive. As we have described in the Introduction, GWs are a dynamical 
effect, while the gravitational redshift is a static effect. Attempting to understand the action of GWs on 
matter with the gravitational redshift is like trying to understand electromagnetic waves using only the 
scalar potential. Next, the fundamental difference between geodesic motion and geodesic deviation motion 
is not appreciated. In stationary spacetimes, such as that surrounding the Earth, a global coordinate system 
can be constructed, and the motion of test particles can be described using the geodesic equation ^^I- In 
non-stationary spacetimes, such as when a GW is present, such a coordinate system cannot be constructed, 
and only the relative motion of test particles can be described, resulting, in the long-wavelength limit, in 
the geodesic deviation equation of motion. 

Among the LIGO community, on the other hand, it is often thought that the use of atom interferometry 
to detect GWs simply involves replacing the light beams in LIGO with atomic beams. The gain in sensitivity 
is then only due to the shorter de Broglie wavelength of the atoms compared with the wavelength of the 
light currently used LIGO; it is still the motion of the mirrors that will be measured. This viewpoint also is 
too naive. Atoms are not photons, and their response to GWs is different than that of photons. As we shall 
see below, the nonrelativistic atoms used in MIGO are moving so slowly that it is the effect of the GWs on 
the speed of the atoms that will be measured; the motion of the mirrors are a secondary effect. 

We hope that the pedagogical nature of this paper will dispel these misconceptions. 

3 Background and Review of Research 

3.1 Matter-wave Interferometry 

We begin with a brief review of the physics underlying matter-wave interferometers. Detailed reviews can 
be found in jIHlEj . 

3.2 Principles of Matter-wave Interferometry 

Matter-wave interferometry is based on the particle-wave duality principle of quantum mechanics, which 
states that every massive quantum object characteristics of both a particle and a wave. Which of these two 
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characteristics it exhibits depends on the properties an experimentalist wishes to measure or exploit. Using 
this principle, various interferometers have been constructed with a wide variety of sources, beam-splitters, 
mirrors and detectors (see ^^lESl)- Following the COW [5] experiment in the 1970s, in the 1990s, Pritchard's 
demonstration of atom interferometry |27j using sodium atoms was followed shortly thereafter by the work 
of Chu and Kasevich |46j with cesium atoms, to perform extremely sensitive measurements of the local 
acceleration due to Earth's gravity. Subsequently, the phase shift of atoms caused by the Earth's rotation 
was measured not only with traditional Mach-Zehnder-type interferometers [3311211, but also with *He and 
■^He superfluid Josephson- Anderson junctions |^ 0SI as well. More recently, matter-wave interference of 
fuUerene (Ceo) molecules has also been reported 

Because of the particle- wave duality, every particle possesses a deBroglie wavelength Ads, and, due to its 
wavelike nature, massive particles can be diffracted, reflected, and coherently beam-split — all the properties 
of light that are needed for light-wave-based interferometry. However, unlike light, A^s depends explicitly 
on the speed of the particle, and its wavelength can be altered continuously. Just as important, a slowly 
moving massive particle responds to the presence of gravity much more sensitively than light, as we shall 
see below. 

Like light-wave interferometry, matter-wave interferometers can be divided into three distinctive parts: 
the source emitting the interfering particle, the 'atom optics' consisting of beam-splitters and mirrors, and 
the detector. In a generic matter-wave interferometer a source emits particles — either in a continuous 
stream or in bunches — which is then split into two (or more) paths at the beam-splitters. (It is important 
that the particle beam is split coherently so that it is not possible — due to the superposition principle of 
quantum mechanics — to determine which path the particle will travel along.) Mirrors are then used to 
change the direction of the beams so that they are eventually recombined and detected. Like the beams 
of light in light-wave interferometers, the particles passing through the interferometer picks up a phase 
shift A(/) = (5j/ — S2l) jh, but now 5^/ and are the action for a particle travelling along the two 
quasi- classical paths 71 and 72 of the two arms of the interferometer. 

3.3 Supersonic Sources 

Supersonic sources were first developed in the 1960s by chemists for the study of chemical reactions (see 
Chapter 2 of 03), and they have a number of properties superior to the effusive sources (ovens) that 
are more often used by physicists. While the velocity distribution of atom beams from effusive sources is 
essentially a Maxwellian one that is fixed by the temperature of the oven, the velocity spread of supersonic 
sources are much narrower. Fractional velocity spreads of 1% — 0.3% for helium — has been achieved from 
continuous supersonic sources, and even smaller fractional velocity spreads can be achieved using pulsed 
sources. Nearly monocnergetic beams with very high intensities can thus be formed. Pritchard, for example, 
has produced 10^^ atoms/cm^ /s/sr sodium beams |271, and Toennies has produced very cold helium beams 
with 1.5 X 10^^ to 1.5 X 10^*^ atom/sr/s [HOI. Throughputs of atoms, measured in atoms per second, in atom 
interferometers using supersonic sources are thus orders of magnitude larger than those obtainable using 
either effusive sources, such as those used by Kasevich |2S1EE1) or magnetic-optical-trap sources, such as 
those used by Chu [45j . In addition, by seeding the beam with the appropriate type of atom, the centreline 
velocity of the atomic beam from a supersonic source can be increased or decreased. This flexibility is 
another advantage that supersonic sources have over effusive sources. 

3.3.1 Continuous Supersonic Sources and Optical Molasses Collimation 

A typical continuous supersonic source functions as follows: A jet of gas from a high-pressure reservoir 
escapes supersonically in free expansion through a nozzle, consisting of a small orifice typically 10 to 100 
microns in diameter, into a differentially pumped low-pressure chamber that has a larger orifice at its output 
end called the 'skimmer'. This skimmer has the appropriate geometry so that it can skim away the hotter 
outer components of the rapidly expanding gas jet, thus leaving only the intense, low-temperature, central 
component of the atomic beam to enter into another differentially pumped chamber. Importantly, after 
expansion, the collision times of atoms in the beam reduce dramatically, and effectively they no longer 
collide with one another. The beam is often further collimated by a slit at its output end, and it could be 
coUimated yet again using a second slit at the output end of yet another differentially pumped chamber 
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before it enters the main vacuum chamber containing the atom-based device, such as an interferometer 
(see for example [HI). With successive stages of differential pumping, by means of either diffusion or 
turbomolecular pumps, one can maintain an ultra-high vacuum in the main chamber that is often needed. 
In addition, using optical molasses techniques |51| . beam collimation can be achieved without throwing 
away most of the atomic beam, as is currently done when collimating slits are used. The combination of 
2D optical molasses collimation with effusive sources has already been achieved by |^ 12^1 ■ Its combination 
with supersonic sources, as we propose here, can yield the high-brightness atomic beams necessary for a 
good signal-to-noise ratio in MIGO. 

3.3.2 Pulsed Supersonic Sources and Optical Molasses Collimation 

Most pulsed supersonic sources function in much of the same way as continuous supersonic sources, but 
with the addition of a fast-acting valve to pulse the beam. Duty cycles of roughly 10"'* can be achieved, and 
as a result, very high fluences can be obtained. Because of the short duty cycles, smaller vacuum pumps 
than those used in continuous supersonic sources can be used instead. Of particular interest is the pulsed 
source developed by Powers et. al. [221 1 which uses laser ablation of metallic sources such as lithium (see 
|53j for a review of laser ablation) to generate the atomic beam. The high-pressure gas chamber is replaced 
by the laser plus metal assembly. As before, optical molasses can be added to reduce the velocity spread, 
and thereby coUimate the beam. 

3.4 LIGO and the Detection of GWs 

With the advantage of being a scalable design, the great majority of the current experimental searches for 
GWs are based on laser interferometers. A number of research groups located throughout the world jSJ] 
are expecting to begin to collect data soon: GEO600, a German-British collaboration; VIRGO, a French- 
Italian collaboration; TAMA300, a Japanese effort; and ACIGA, an Australian effort. LIGO, the US-based, 
international collaboration, is currently collecting data, and LISA, a space-based laser interferometer system, 
is currently in the initial planning stage. We shall focus primarily on LIGO in this paper. 

LIGO is a set of three interferometers based at two locations — Hanford, Washington and Livingston, 
Louisiana — separated by 3020 km [HHl- All three instruments are based on Michelson interferometers with 
Fabry-Perot arms. The physical length of the each arm of the main LIGO interferometer is 4 km, and the 
Fabry-Perot interferometer increases the optical path of the arm a factor of 2B = 150. In its current LIGO I 
configuration, a 8 W Nd:YAG laser is used with = 1.064 /xm. At the end of each arm is a massive mirror 
suspended vertically within a vacuum chamber; the location of the mirror assembly for the end mirror must 
be held in position within 10"^" to 10~^^ m with respect to the central beam splitter of the interferometer. 
As we shall see in the next section, the passage of a GW through the interferometer shifts the position of the 
end mirrors, thereby causing a net phase shift of the light, which is then measured. The laser interferometer 
used by LIGO thus provides the means of accurately measuring the position of the mirrors over a large 
distance. 

While each end-mirror could be fixed to the mirror assembly, the end mirrors would, through their 
connection with the mirror assemblies, be mounted to the ground below. The response of the interferometer 
to the passage of the GW would then depend on the material properties of both the ground and the frame 
of the interferometer across the 4 km spanned by LIGO's arms. Thus, by fixing the mirrors to the mirror 
assemblies, one ends up measuring the response of the Earth to the GW, and not of the mirrors. If this is 
done, sources of systematic errors — such as changes to the material properties in the Earth along LIGO's 
arms — cannot be controlled. Consequently, instead of fixing the end mirrors to the mirror assembly, they 
are suspended on piano wires, and are thus decoupled as much as possible from sources of uncontrollable 
systematic errors. Even so, at frequencies below 125 Hz thermal noise in the piano wires will begin limiting 
LIGO's sensitivity, and below 40 Hz, seismic noise causes a rapid decrease in sensitivity. At frequencies 
above 125 Hz, shot noise begins limiting the phase sensitivity of the interferometer. 

Construction of LIGO began in 1996, and the main interferometers were commissioned in 2001 [SH]- The 
first science runs were started in June of 2002, and these data are currently being analysed. The installation 
of LIGO II — designed to be used for GW astronomy — is expected to begin in 2006. 
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Figure 5: Schematic diagram of the horizontal MIGO configuration, with a GW incident normal to the 
plane of the interferometer. Only the x polarization contributes to the phase shift; the + polarization does 
not. Diffraction orders other than n = ±1 are left out for clarity. 

4 A Matter-wave /nterferometric Gravitational-wave 
Observatory (MIGO) 

In this section we outline two possible configurations of MIGO and the technologies that could be used to 
construct them. We also derive the phase-shifts expected for this MIGO, and use this calculation as a guide 
in determining the specifications of MIGO based on current technology. 

4.1 Theoretical Basis 

In this subsection we outline the physics underlying LIGO and MIGO, and demonstrate that they are both 
based on the same underlying physics, one at the classical level, the other at the quantum level. The 
treatment here will be semi-quantitative, and will focus on the fundamental physics. A complete derivation 
of MIGO's phase shift is given in Appendix B, and the reader is referred to 7 for the derivation of LIGO's 
phase shift. Like LIGO, we construct our coordinate system along the lines given in the Introduction (see 
|12| as well), and we measure the position and movements of particles with respect to the initial beam 
splitter (see figures |5l and Ell . 

The basic physics underlying LIGO is as follows: Since the mirrors at the end of each of LIGO's arms 
are hung as pendula with piano wire, they are free to oscillate along the interferometer's beam line. As a 
GW passes through the interferometer it slightly shifts the position of both mirrors by a small amount Ax. 
This in turn causes a shift Ac/jljqo ~ ^■nBAx/X-f in the phase of light (with wavelength A^) used in the 
interferometer [7||SH]. (The restorative force of the pendula on the mirrors can be neglected here.) It is thus 
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Figure 6: Schematic diagram of the vertical MIGO configuration, with a GW incident normal to the plane 
of the interferometer. Only the x polarization contributes to the phase shift; the + polarization does not. 
Diffraction orders other than n = ±1 are left out for clarity. 

the shifts in the position of the mirrors caused by the GW that is being measured, and not the changes 
caused by the GW in the properties of the light, as it is often thought among the AMO community. The 
light used is only a very accurate means of measuring shifts in the mirror's position. 

As described in the Introduction, the motion of the mirrors in the interferometer is govern by the geodesic 
deviation equation Eq. . Since the passage of the GW causes only a small change in the position of the 
mirror, Aa; « LLicoh+fi where Lligo = 4 km is the physical length of an arm of LIGO, and h+ is the 
amplitude of the + polarization of the GW; the x polarization does not affect the motion of the mirrors 
appreciably. This shift in the mirror's position changes the optical path length of the interferometer, and 
causes a phase shift in the light used in the interferometer given by AcpLiao = 271 BLLiGoh+/^-y [SEIIZI- 
The factor B is the number of round trips that the light beam makes inside the Fabry-Perot interferometers 
in LIGO's arms. Strictly speaking this expression for AcpLiGO is only applicable for GWs with frequencies 
~ 125 Hz; at 4 km, the arms of LIGO is equal to the reduced wavelength of a GW in the upper end of its 
frequency response spectrum, and causality has to be taken into account (see |S7], page 550). While design 
characteristics of LIGO, such as power recycling and shot-noise modifications by the Fabry Perot arms, will 
change the sensitivity from this simple form, it is accurate enough for our purposes in this section. 

If the phase shift for LIGO is proportional to the length of its arms, how, then, does the phase shift 
for MIGO depend on its size? Figures [S] and are drawings of two possible configurations for MIGO with 
the force-field lines of a x polarization GW drawn in at the origin of our coordinate system. A complete 
description of these configurations will be given in the next subsection. For now, we note that on a classical 
level, the motions of the atoms used in MIGO are governed by Eq. Q as well. However, while in LIGO 



14 



one considers shifts of a mirror's position placed at a distance from the central beam splitter (which is 
the origin of LIGO's coordinate system) due to the passage of a GW, in MIGO one considers shifts in an 
atom's velocity as it travels through the interferometer due to the GW. Thus, while an atom may have an 
initial velocity Vq after it passes through the first beam splitter (in contrast with the end mirror in LIGO 
which has an initial position), the passage of a GW will cause this velocity to be shifted by an amount 
roughly Aw* ~ Wp/iy r/2 after it has traversed the length of the interferometer in a time T. With a typical 
\hij\ ~ 10"^*^, this velocity shift of the atoms is extremely small, and must be measured on the quantum 
level by means of matter- wave interferometry. 

Note that the tidal force-field lines sketched in figures |S1 and IHl break the bilateral symmetry of the 
interferometers. It is this breaking of the bilateral symmetry by a x polarization GW that leads to a 
non-zero phase shift. 

Changes in the atom's speed \Av\ result in changes to its deBroglie wavelength, and hence to cumulative 
changes in its quantum phase. Indeed, since the change in the atom's position due to this speed shift 

after it passes through the interferometer is roughly ALmigo — 2|Aw*|r ~ LmigqT hij (the factor of 2 
comes from having two arms in the interferometer), where Lmigo is the effective size of MIGO, and the 
transit time T Lmigo/\vq\- Like a laser interferometer this change in position will cause a phase shift, 
A4>MiGO = '^T^ ALmigo / ^dB ~ "^^m/goI^u 1/^' where the deBroglie wavelength of the atom is now used 
instead of the wavelength of light. Thus, in contrast to LIGO's phase shift A(f)]^iQO: which is proportional 
to hij, MIGO's phase shift AcpMiGO is proportional to the rate of change of the GW amplitude hij ~ //ly. 

Although simplistic, this rough derivation of the MIGO phase shift nevertheless elucidates the underlying 
physics of MIGO. The above expression for the MIGO phase shift is close to the exact equation calculated 
in Appendix B after a detailed analysis of the MIGO configurations shown in figures and El 
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|A0A//Go(/)| = 27rjAfh^if)\FifT)\, (5) 

were A — L^uqo is the effective area of the interferometer, m is the mass of the atom, / is the frequency of 
the GW, and is the amplitude of the x polarization of a GW. Equation Q follows from Eq. (jSJ, the above 
expression for A^ljco, taking |i^(/r)| w 1/2 for freely suspended mirrors, and setting the shot-noise limits 
of the phase shift of MIGO equal to that of LIGO. The resonance function F{fT) depends on the specific 
configuration of the interferometer, and measures the resonances between the GW and the interferometer. 
Precise forms for A and F{fT) are given in Appendix B. The corresponding shot-noise- limited sensitivity 
for MIGO is 

l^irsMIGO ^ ^ /gx 

2T:mAf\F{fT)\m/^' ^ ' 

where N is the number of atoms passing though the interferometer per second. In comparison 

-.LIGO _ / ^7 1 ^'"^ 27r/ 



where /q is the power of the laser, and 77 is the photodetector efficiency (Eq. 123a of |7j). Thus, while 
^(f)sh^'^ decreases at higher frequencies, ^(/)^o?'^ increases. This complementarity between the two 
sensitivities is due to a fundamental difference in the signals being measured by LIGO and MIGO: LIGO 
measures the position, while MIGO measures the velocity of test masses. 

Unlike LIGO, MIGO is only sensitive to the x polarization, not the -I- polarization. Also, while the 
phase shift of LIGO scales with its length, the phase shift of MIGO scales with its area; the larger the area, 
the smaller the amplitude of the GW that can be detected. When the transit time is much greater than the 
period of the GW, \F\ ^1/2 for freely suspended mirrors; \F\ ^ nfT when it is much less than the period. 
What is not included in Eq. ^ is the Sagnac effect caused by the rotation of the Earth |2SI , and the 
phase shift caused by stationary sources of curvature such as the Earth. While for the configurations shown 
in figures [S] and these phase shifts are expected to be very large, they are, however, steady-state, and thus 
can be isolated from the time-varying signal caused by a GW. 
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4.1.1 Slow Atoms and the Response of MIGO 

With the characteristic speed of an atom from a supersonic source being in the hundreds of meters per 
second range, and T in the range of seconds, we would naively expect MIGO to respond very slowly to 
the passage of a GW, which moves through the interferometer at the speed of light. (However, even for a 
laser-interferometer-based observatory such LIGO, which has arms that are roughly equal to the reduced 
wavelength of a 10'' Hz GW, the transit time of light through LIGO is comparable to the period of the GW 
as well.) Because of the cumulative effect of the GW on the phase of the atom, we would expect A(/)m/go 
to be the result of averaging the phase of the atom over many periods of the GW, and thus to be close to 
zero. This does not happen for the following two reasons. 

Firstly, half-way through the interferometer the atoms in the beam hit, and are reflected off, a set of 
mirrors. These mirrors impart on the atoms an instantaneous force whose strength is proportional to the 
instantaneous velocity — including the changes to the atom's velocity caused by a GW as it travelled between 
the beam splitter and mirror — of the atom normal to the mirror. From the force-field lines drawn in figures 
El and El the magnitude of this force will be different for atoms travelling along the n = +1 and n = — 1 
paths, and they will induce a net phase shift between the atom travelling along the two different paths. 
This can be seen explicitly in the derivation of the MIGO phase shift given in the Appendix B, where the 
force that the mirrors exert on the atoms results in a jump condition for the atom's velocity in the direction 
perpendicular to the mirror's surface. For high-frequency GWs this effect of the mirror's impulsive force is 
the dominant contribution to A^m/go, and is the reason why \F\ ~ 1/2 for freely suspended mirrors when 
T is much longer than the period of the GW. At lower frequencies, when T is much shorter than the period 
of the GW, the mirror's impulsive effect is very much reduced. 

Secondly, notice from Eq. (O that due to the tidal nature of the GW, its effect on the motion of the atom 
increases as the atom moves away from the beam splitter; the atom 'sees' a larger effective acceleration 
later in its path through than interferometer than it did at the beginning. Although for low frequency GWs 
this results in a relatively small resonance function \F\ ^ T^fT, the phase shift nonetheless does not average 
out to be zero as quickly as one might expect. 

The relative sizes of MIGO compared with LIGO given by Eq. would also seem to be counterintuitive. 
MIGO makes use of slowly-moving, nonrelativistic atoms to make its measurements, while LIGO would seem 
to make use of photons moving at the speed of light c. At first glance it would seem that MIGO should 
be less sensitive to GWs than LIGO by some power of Vatom/c. This, however, would be an erroneous 
argument. As outlined in 3.4, it is not the effect of GW on the light used in the laser interferometer that is 
being measured in LIGO; it is the effect of the GW on the test masses (the end mirrors), which are at rest, 
that is being measured. Indeed, it is precisely because the atoms are slowly moving, i.e., nonrelativistic, 
that the effect of GWs on their motions are much more readily measurable than their effect on light. A GW 
can readily change the speed of nonrelativistic atoms, but it cannot change the speed c of light. 

As a consequence, we would expect the slower the atom is, the more sensitive MIGO will be, as long 
as the atoms in the beam do not decohere because of long transit times through the interferometer. As 
counterintuitive as this conclusion may be, it can be arrived at through the following thought experiment. 
Suppose we replace the slowly moving, nonrelativistic atoms in MIGO by a series of faster and faster atoms 
until the atoms in the beam approach the speed of light. In this ultrarelativistic limit, the atoms behave 
much like the light used in LIGO, and it is apparent the horizontal MIGO configuration would be similar 
to LIGO, but with mirrors that are rigidly attached to the frame of the interferometer, instead of the freely 
moving mirrors in LIGO. Since it is well known that the sensitivity of LIGO decreases greatly if rigidly 
mounted mirrors are used instead of suspended ones, MIGO's sensitivity to GWs when ultrarelativistic 
atoms are used is expected to be much worse than if slowly moving atoms are used. Seen another way, 
because the speed of ultrarelativistic atoms are so close to the speed of light, changing their speed by even 
a small amount requires a great deal of energy. The speed of nonrelativistic atoms, on the other hand, 
are much more readily changed, and for the same amount of energy, a GW can shift the speed of a slow, 
nonrelativistic atom much more than the speed near c of an ultrarelativistic one, and thereby can cause a 
much larger phase shift for the nonrelativistic atom. To summarize, we would expect that the sensitivity 
to detecting GWs to decrease if ultrarelativistic particles — which behave like photons — are used in MIGO 
instead of slowly moving, nonrelativistic atoms. However, the slowly moving atoms must not decohere due 
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to collisions during the transit times through the interferometer. 



4.1.2 Comparing MIGO with Other Approaches to the Detection of GWs using Matter-wave 
Interferometers 

Equations Q), ||SJ|, and the conclusion that the sensitivity of atom-based interferometers to GWs is larger 
in comparison to light-based interferometers, are also surprising when compared to the results of Linet- 
Tourrenc ^28^ , Stodolsky , Cai-Papini , and Borde and coworkers 35i |36j . Complete criticisms of 
these approaches will be given elsewhere We note here simply that the approaches taken by both 

Linet-Tourrenc and Stodolsky have their roots in the geodesic equation of motion. Linet-Tourrenc used 
the quasi-classical approach based on the relativistic single-particle Hamiltonian for a relativistic particle 
derived from the geodesic equation; Stodolsky based his approach on the relativistic single-particle action. 
Both authors do not explicitly define their coordinate systems. They both find a phase shift of the form 

f dt 

^4>geodesic ^ ~^ J ^vP^'P'^' (8) 

from the equation above Eq. (6.1) of and Eq. [3.2.3] of j2S- (To conform with standard units, we have 
restored the requisite factors of c and h.) The integral is over a closed path 7 in spacetime that the particle 
takes through the interferometer, pi is the momentum of the particle, and E is its total energy including 
the rest mass. In comparison, if we also, for the moment, neglect the action of the GW on all parts of the 
interferometer (which is an unphysical assumption), we find the change in the action to be 

m /" / 1 ^, 1 dhij ■ ■\ 

^(PgeodeviaUon ^ ^ * (^2^" ^ 2~dt^ ) ' 

The difference between ^(pgeodesic and A.(j)geodeviaUon is readily apparent, especially when the nonrelativistic 
limit of geodesic is taken. ^4>geodesic depends on the momentum of the particle and is independent of 
the length of the path travelled. Because of this momentum dependence of l^^cj) geodesic, it is expected that 
matter-wave interferometers will become more sensitive the closer to the speed of light the atom moves; 
i.e., the more ultrarelativistic the particle becomes, the larger its phase shift. This can be seen explicitly 
when the nonrelativistic limit of Stodosky's expression (Eq. (6.1)) for the phase shift of his interferometer is 
taken, in contrast to our expression for the phase shift of MIGO. On the other hand, A(j)geodeviation increases 
linearly with the separation a;*, as it should, and the larger the interferometer, the larger the change in the 
phase, in agreement with the tidal nature of GWs. This linear scaling of the phase shift with the size of 
the interferometer is also a result on which LIGO is based. In addition, unlike previous studies, the effects 
of the mirrors on the atom's phase shift was also considered in our analysis. As a consequence, the phase 
shift of MIGO calculated in this paper has a dependence on the velocity of the atoms in the interferometer 
that is opposite from what Stodolsky calculated, when the nonrelativistic limit of his Eq. (6.1) is taken. 

It should be emphasized that the differences between A(j)geodesic and A(j>geodeviation is due to more than 
the fact that one was derived in a nonrelativistic framework, and the other within a relativistic one. This 
can be seen by comparing the equations of motion for the particle used in the two approaches. Both Linet- 
Tourrenc and Stodolsky's approaches use the geodesic equation of motion, which, when GWs are present, 
becomes 

d^r ] dh ■ 

dt^ ^ 2 dt ' ^ ' 

instead of the geodesic deviation equation Eq. It is, however, well known that GWs affect all parts of 
the system, and since only differences in positions can be measured, it is the geodesic deviation equation 
that should be used, and not the geodesic equation. This becomes apparent in Eq. H1Q(I with an acceleration 
for the particle dependent on the derivative of /lij, and thus, from Eq. H15(l . on the connection. Equation 
(|1(J|I is thus a frame- dependent quantity, and will only have meaning only when a coordinate system is 
explicitly chosen. This choice was not explicitly made either by Linet-Tourrenc or by Stodolsky |58j, nor 
by Cai-Papini or by Borde and coworkers. This coordinate choice is explicitly made, and is well defined, in 
the derivation of the geodesic deviation equation of motion, Eq. Q, and is the reason why the acceleration 
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of the particle here depends on the Riemann curvature tensor, a bame-independent object. These are the 
underlying physical reasons why the above authors obtained the incorrect expression for the phase shift. 

On the quantum-field-theoretic level, Borde and coworkers start with the Dirac equation in curved 
spacetime. They used the standard formalism for quantum field theory in the presence of linearized gravity 
to study the phase shift that a GW induces on a generic atom interferometer using nonrelativistic, spin-half 
atoms. They obtained an expression (Eq. (92) of (HSI) for the quantum phase shift of the atom that has 
a power-law dependency on X±dB and Xgw that is very much different than ours, and instead is similar 
to Linet-Tourrenc's and Stodolsky's. This also led them to the erroneous conclusion that atom-based 
interferometers are at best no more sensitive to GWs than light-based interferometers. Recall from Eq. (|3J), 
that the tidal effects of a GW on an atom increases with the separation distance between the observer and 
the atom. Borde's governing Hamiltonian (Eq. 89 of [321), on which their analysis was based, does not 
include these tidal effects of GWs, resulting in an incorrect expression for the quantum phase shift for the 
atom. 

4.2 Description of MIGO Configurations 

The two configurations of MIGO shown in figures Eland El represent two possible operating extremes. In the 
horizontal configuration shown in figure O atoms are emitted from a continuous supersonic source with a 
velocity high enough that the acceleration due to Earth's gravity does not appreciably alter an stom's path 
by the time it traverses the interferometer. Indeed, not shown in the figure is a Zeeman laser accelerator^ 
working oppositely to a Zeeman slower [SI], that can accelerate the alkali atoms to even higher velocities. 
The atomic beam from the supersonic source is automatically collimated, and the distribution of velocities 
of the atoms in the beam transverse to the beam typically has a Avt/vs = 1% for most atoms, and can 
be decreased to 0.3% for helium j45| . Alkali atomic beams from supersonic sources can then be further 
collimated by means of 2D optical molasses. The atom beam then passes through the first beam splitter, a 
nanofabricated transmission diffraction grating which, as with the diffraction of light, splits the beam into 
different diffraction orders n. 

As shown in figure|3 only the n — ±1 orders are used here. In is important to note that the diffraction 
grating splits the atomic beam coherently. Consequently, from the superposition principle of quantum 
mechanics, it is not possible to determine which path any one atom in the beam will take; there is a finite 
probability amplitude it will take the n — +1 path, and due to the bilateral symmetry of the interferometer 
along the horizontal line, the same probability amplitude that it will take the n = —I path. Irrespective of 
which path is taken, the most probable trajectory between the initial beam splitter and the mirror for the 
atom is determined by Eq. J^J (see j43| for the use of the Feynman path integral in atom interferometry). 
The mirror then refiects the atom by exerting an impulsive force whose strength is proportional to the 
velocity of the atom normally incident upon the mirror. (See |59) for recent advances in atom optics, such 
as coherent atomic mirrors made from crystals, and the possibility of using curved mirrors to coherently 
focus and coUimate atomic beams.) Once again the most probable trajectory of the atom between the mirror 
and the final beam splitter — where the n = +1 path and the n = —1 path is recombined — is determined by 
Eq. 0; the atoms are then detected using standard methods jlHl- In the absence of gravitational effects, 
the bilateral symmetry of the interferometer ensures that the n = -1-1 path will be the same as the n = —1 
path, and no phase shift will be measured. This symmetry is broken when a GW is present, however; from 
the force-field lines drawn in figure [3 we see that path that the atom takes for the n — +1 order will be 
slightly different than the one taken by the n = —1 order, leading to MIGO's phase shift. 

For the horizontal MIGO configuration, A = L±L\^/2 — the actual area enclosed by the interferometer — in 
Eq. (jSJ (see Appendix B). Since the atoms travel in approximately straight lines through the interferometer, 
the length i|| of this MIGO configuration is determined by its width L±, and the ratio of the horizontal 
velocity of the atom to its transverse velocity v± after the initial beam splitter (see Appendix B). The 
sensitivity of a horizontal MIGO is thus proportional to i^. 

In the vertical configuration for MIGO shown in figure 13 alkali atoms are emitted from a pulsed super- 
sonic source, and are then slowed down using a Zeeman slower |51| to velocities so slow that the acceleration 
due to Earth's gravity g now dominates the trajectories of the atoms through the interferometer, similar to 
the atom fountain geometry of Chu and Kasevich (46|. The beam is collimated using 2D optical molasses 
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as well 1201 ■ The velocity of the atoms Vs along the beam is slowed to such an extent that after passing 
through the initial beam splitter the atoms in the n = ±1 orders now traverse an almost parabolic trajectory 
to mirrors placed at the maximum height of their trajectories. The atoms are once again reflected at the 
mirrors, and fall back downward in an almost-parabolic trajectory where the paths are recombined by the 
same beam splitter, and subsequently detected. Once again, the bilateral symmetry of this interferometer 
along now the vertical line ensures that the gross motion of the atoms on the left side of the interferometer 
is reflected on the right side. The usual parabolic trajectories of the atoms is shifted slightly by an amount 
determined once again by Eq. Q due to the passage of a GW, and will again be different for atoms travelling 
along the left hand path instead of the right hand path. It is this asymmetrical shift by the GW in the 
atoms trajectories that leads to the overall phase shift. 

Similar to the horizontal MIGO configuration, A = L±L\\ for the vertical MIGO configuration as well 
(see Appendix B), but because the acceleration due to gravity g slows the atom down, it spends more time 
in the presence of a GW inside the interferometer than the atoms in the horizontal MIGO configuration do. 
Consequently, L|| = gL\/Sv\, where the horizontal velocity of atom after the beam splitter, is inversely 
proportional to the periodicity of the diffraction grating used. Thus, the area A is proportional to the cube 
of the width of the interferometer, and the sensitivity of the vertical MIGO configuration is proportional to 
L\, and not to L\ as for the horizontal MIGO configuration, or even to L± as it is for LIGO. 

Like LIGO, both MIGO configurations make use of mirrors to redirect the atomic beams to beam 
splitters, and like LIGO, the response of MIGO to GWs will, to a certain extent, depend on how the mirrors 
are attached to the frame of the interferometer, and thus how the interferometer is attached to the Earth. 
However, because it is effects of the GW on the interfering particle, i.e., the atoms, that is being measured 
in MIGO, and not the mirrors as it is in LIGO, MIGO is much less dependent on how this attachment 
is made than LIGO is. Indeed, from Appendix B, we found that for fT >> 1, the response function 
1/2 < F{fT) < 1, irrespective of how the mirrors are attached to the interferometer. In fact, the lower 
limit of 1/2 that is reached if the mirrors are freely suspended as in the case of LIGO, and the upper limit 
of 1 is reached if the mirror was firmly attached to an interferometer frame that can be approximated as 
infinitely rigid [Snj. MIGO thus becomes more sensitive the more firmly the mirrors are attached to the 
interferometer frame. This is in direct contrast to LIGO, which become less sensitive. Indeed, if like MIGO 
we model the connection of LIGO's mirrors to the interferometer frame as a spring with resonance frequency 
/o and loss factor Q, we find from Eq. (37.16) of 0] that the shift in position of the end mirrors will be 

h+L, (11) 



P - fi + *//o/Q 

In the limit where /o >> /, when the end mirrors can be approximated as mounted to an interferometer 
frame that is infinitely rigid. Ax « 0, and the LIGO looses all sensitivity to GWs. 

From the sketch in figure |H| we see that the vertical MIGO configuration is a combination of Chu- 
Kasevich's atomic fountain interferometer with a Pritchard-type interferometer. Previous to Chu and Kase- 
vich's work, Zacharias (see |HJ|) attempted to construct an atomic fountain in the 1950s and failed. Because 
of the effusive source used, there was a broad distribution of velocities in the atomic beam, and as the atoms 
in the beam climbed up the gravitational potential of the Earth, they started to slow down, and other, 
faster atoms started to overtake them. Scattering with the faster-moving atoms in the beam took place 
and coherence was lost. Chu and Kasevich avoided this by using laser-cooled atomic beams with narrower 
velocity distributions, as does MIGO. 



4.3 Potential Sensitivity 

The effectiveness of both MIGO configurations in detecting GWs can be seen in figure [7| where the char- 
acteristic amplitude and frequencies of GWs emitted from burst sources (the strongest emitters of GWs) 
such as supernova explosions and the coalescence of black holes, are plotted along with the amplitude and 
frequencies of the smallest of these types of signals that are detectable by MIGO, LIGO (both the current 
LIGO I configuration and the LIGO II configuration slated for operation in 2007), and LISA. Two of the 
configurations — the Earth-based vertical and horizontal — was designed to have a sensitivity equal to the best 
sensitivity of LIGO in its frequency range, and the third configuration — the Space-based horizontal — to have 
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a sensitivity equal to the best sensitivity of LISA in its frequency range. The signals from various classes of 
burst signals were replotted from Thome's figure 9.4 jTj; the specifics of these sources can be found therein. 
The plots LIGO I, LIGO II and LISA have been updated from Thome's original figure, and were generated 
from [7], inn J and [53], respectively, using Eq. Ill of \7\. 

All three graphs for MIGO were generated using Eq. (|SJ), and a. N = 10^* atoms per second (as compared 
with 8 X 10^^ photons per second for LIGO). While N = 10^^ atoms per second is very large, it is important 
to note that what is important is the number of atoms per second through the interferometer, and not the 
number density. Indeed, if the number density of atoms is too high, the atoms in the beam will condense, 
and beam coherence will be lost. By suitably constructing the supersonic source, and by choosing the proper 
operating regime for it, it may be possible ensure that the 10^^ atoms per second needed to achieve figure[7| 
is spread over a large enough cross section of the beam to prevent condensation. We will comment further 
on this later in 4.4. 

A nanofabricated transmission diffraction grating with a periodicity of 32 nm was used in calculating 
figured Although only reflection diffraction gratings with this periodicity have been made by the Nanoelec- 
tronic Research Centre at the University of Glasgow by flipping the source and the detector array in 
figures Inland El the MIGO configurations can be modified to use refiection gratings if transmission gratings 
cannot be constructed. As with LIGO, because GWs cannot be shielded, the passage of a GW through 
MIGO will shift all its parts, including the mirrors, and the response of the mirrors to the GW introduces 
an additional shift in the velocities of the atoms. However, since the mirrors in MIGOs do not have to be 
free to move, they can be constructed so that the effect of the GW on their motion is critically damped, as 
they were taken to be when figure was graphed. 

For the Earth-based vertical MIGO, bosonic ^'^^Cs atoms from a 2000 K continuous supersonic source 
are used in an interferometer with a width Lj_ = 1.3 m and a height of L|| — 237 m. Zeeman cooling of 
the vertical velocity of the atoms from the beam is used to decrease the velocities of the atom by a factor 
of 11.6, and additional 2D optical molasses cooling is used to narrow by a factor of lO'' the velocity of the 
atoms perpendicular to the beam. The Space-based MIGO, which is depicted as a horizontal configuration, 
also uses ^'^^Cs atoms emitted from a 2000 K continuous supersonic source, but with 2D optical molasses 
cooling used only to collimate the beam. As a result, while L± = 2 m, the length of the interferometer is 
L|| = 17 km, as compared to 5 million km for LISA. (Although a throughput of 10^^ atoms/s is high for a 
space-based MIGO, atoms can be recycled from detector to source.) The final Earth-based horizontal MIGO 
graph also has a L± = 1.2 m, but uses fermionic ^Li atoms from a 2000 K supersonic source in conjunction 
with a Zeeman laser accelerator to increase the beam- velocity of the atom by a factor of 1.86. By angling 
the beam slightly upward, the maximum height of the beam as it travels across the length i|| = 4 km of 
the interferometer is only 0.41 m; 2D optical molasses is used to collimate this beam as well. This MIGO 
configuration, while not optimal, is designed so that whole interferometer fits within one of the 1.2 m wide 
evacuated beam tubes that form the arms of LIGO. Enlarging sections of the beam tube to ^ 2 m will 
allow ^^'^Cs to be used instead of ^Li, resulting in a 37-fold increase in the Earth-based horizontal MIGO 
sensitivity as shown in figured 

In figuredwe see that by replacing the lasers in one of LIGO's arms with the ^Li horizontal MIGO, the 
sensitivity of the GW observatory can be increased by a factor of 10. At 1.2 m x 237 m, the Earth-based 
vertical MIGO has the same sensitivity as LIGO II in its most sensitive frequency regime, while at 2 m 
X 17 km, the Space-based MIGO has the same sensitivity as LISA in EISA's operating regime. But all 
three MIGO configurations are sensitive to a much wider range of GW frequencies than LIGO and LISA, 
and they extend the range of GW observatories into frequency ranges — in particular between 0.1 and 10 
Hz — not currently accessible by LIGO or LISA. In particular, note the sensitivity to high-frequency GWs 
of all three MIGO configurations graphed in figure d and the gradual insensitivity of LIGO to frequencies 
above 125 Hz and of LISA to frequencies above 2 x 10^^ Hz. This is due to the difference in the signal 
being measured by matter-wave-based and light-wave-based interferometers: MIGO measures shifts in the 
atom's velocity, and thus the rate of change of the GW amplitude with time; the larger the velocity shift, 
the larger the phase shift. This rate of change increases at higher frequencies, and is the underlying reason 
the smallest-amplitude GW that three MIGO configurations can detect in figured (decreases as 1//^^^ at 
high frequencies. In contrast, LIGO (and LISA) measures the position of the end mirrors, which depends 
only on the amplitude of the GW, and only indirectly on its frequency. For the frequency of the GW to be 
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Figure 7: The ability of MIGO to detect GWs from various classes of burst sources is plotted, and compared 
with that of LIGO I 7 , LIGO II 62 and LISA [S3] . All three MIGO configurations are sensitive to a broad 
range of GW frequencies; they thus extend the range of GW observatories into frequency ranges not currently 
accessible. 
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measured, the smallest-amplitude GW that LIGO can detect must increase as /^/^; it actually increases as 
/•^/^ because of shot-noise limits to A0l/go- Indeed, if we compare shot-noise-limited sensitivities, 
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where N-y is the rate at which photons enter LIGO. 



4.4 Design Considerations 

The inherent difference in the signal that they measure — MIGO measures velocity while LIGO measures 
position — impacts their constructions, and the fundamental limits to the size of the phase shifts both are 
capable of measuring. Both LIGO and MIGO will be subject to systematic errors arising from environmental 
perturbations, such as thermal and seismic noise; stray electric and magnetic field gradients; and fluctuations 
in g. Because the temperature of the beam splitters and the mirrors in both MIGO and LIGO cannot be 
zero, temperature-driven changes in their positions and velocities will induce thermal noise into phase shift 
measurements. Random seismic movement of the ground under the interferometer will introduce vibrational 
noise as well. However, while thermal noise starts limiting LIGO's sensitivity below 125 Hz and seismic 
noise effectively cuts it off below 40 Hz, their effect on MIGOs sensitivity is not nearly as great. Thermal 
and seismic noise have such a large effect on LIGO because LIGO's end- mirrors must be as close to a 
freely- falling mirror as possible. MIGO, on the other hand, makes use of beam splitters and mirrors that 
can be solidly mounted on an interferometer frame whose material properties can be well characterized 
and controlled, because of the relatively small size of MIGO. Standard vibrational isolation technologies 
can thus be used to minimize seismic noise to levels not attainable by LIGO; for these reasons seismic 
noise limitations are not shown in figure lO. However, temperature-driven oscillations on the surface of 
the mirrors — even though they arc rigidly mounted — will induce random velocity (Dopplcr) shifts to the 
atoms, and subsequently random phase shifts. These oscillations, although not important for LIGO, will 
slightly affect the sensitivity of MIGO at very low frequencies. The effect of these oscillations on MIGO's 
sensitivity is estimated in Appendix C, and are included in figure [T] Unlike LIGO, the atoms in MIGO 
are subject to stray electric and magnetic fields, however, these effects can be made negligibly small with 
proper shielding and design. Fluctuations in g are due to density fluctuations caused by seismic waves, and 
will limit the low-frequency response of LIGO i22|- Its effect on MIGO is unknown, and can be eliminated 
with a space-based system. 

At a fundamental level, unlike the photons used by LIGO, the atoms used in MIGO can interact among 
themselves; this requires a certain amount of care to be taken in the design of MIGO. Indeed, fcrmionic 
^Li was chosen for the Earth-based horizontal MIGO precisely because as a fermionic atom it will not 
scatter with other ^Li atoms as readily as ^Li, a bosonic atom, will. In this sense, cold, fcrmionic ^Li atoms 
will behave more like the noninteracting photons in LIGO. Moreover, number squeezed states for fermions, 
which are very robust, can greatly reduce the shot noise of the atomic beam [39], [23]. This correspondingly 
reduces the required throughput of the atoms needed for the same sensitivity. 

Any collision of the atoms within the beam — either with other atoms in the beam, or with stray atoms in 
the interferometer — will introduce random noise in Ac/jmigO: ^^^d decrease the contrast of the interference 
fringes. The cross-sectional area of the beam has to be made large enough, and the transverse velocity 
spread of the beam — narrowed using 2D optical molasses — has to be small enough that this does not occur 
within the transit time of the atom through the interferometer. Scattering with the background gas can 
be made negligible by using ultra- high- vacuum systems. Next, the supersonic source must be run at a 
high enough temperature that condensation of the gas of atoms into clusters — either solid or liquid — does 
not occur. At the same time, the pressure of the supersonic source must be high enough to generate a 
high intensity of atoms through the interferometer. Since we are dealing with atoms which can condense 
into liquids and solids, the allowed running temperature and pressure of the source will depend on the 
phase diagram of the atom, and the degree of supercooling of the gas that takes place. These issues are 
well studied for helium 65 (although they are not as well studied for the alkali atoms), and in general 
there does not seem to be a fundamental obstacle to finding a set of viable operating parameters. What is 
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more challenging is the throughput of 10^^ atoms per second used in generating figure [7| which is close to 
the limits of current supersonic-source technology. However, because of the sensitive dependence of MIGO 
sensitivities on its size, a decrease in this throughput by a factor of 100 can be made up for by increasing the 
L±_ of a horizontal MIGO by a factor of lO^/^ = 3.16, or that of a vertical MIGO by a factor of lO^/^ = 2.15; 
decreases in the atom-beam throughput can thus be readily compensated for by modestly increasing the 
size of the interferometer. 

At frequencies below roughly 125 Hz the sensitivity of LIGO decreases due to thermal noise, and radiation 
pressure on the end mirrors j62| . Photons in the laser impacting the mirrors not only exert pressure (and 
thus move) the mirrors, but because they are not perfect, the mirrors will be heated by the incident beam. 
Both factors will cause an overall decrease in LIGO's sensitivity. One would expect that since 10"'^^ atoms 
per second should pass through MIGO to achieve high sensitivity, the impact pressure of the atoms on the 
mirrors, combined with the heating of the mirrors caused by the inelastic scattering of the atoms, would 
limit the sensitivity of MIGO as well. They do not for the follow reasons. 

Quantum interference places stringent limits on the decoherence of the atom beam as it is diffracted 
from the beam splitters, and reflected off the mirrors. (Scattering of atoms in the beam off of background 
gas can be shown to be negligible.) In particular, it requires that there exists no possibility of 'which-path' 
information for the atom inside the interferometer. This in turn requires that diffraction of atomic beams 
off of beam splitters, and their reflection off of mirrors be elastic and coherent. Consequently, an atom can 
only change the centre- of -mass momentum of the beam splitters or mirrors, and cannot deposit energy into 
them. Since the mirrors and beam splitters of MIGO can be firmly mounted to the frame, the effects of this 
centre-of-mass momentum transfer can be minimized. 

A measure of the inelastic versus elastic components of these processes is based on the Debye- Waller 
factor W through the intensity ratio ///q = exp(— 2iy), where / is the diffracted (or reflected) intensity and 
/o is the incident intensity. The factor is a measure of the fluctuations of the atoms in the crystal that 
diffracts the incident atomic beam. For diffracted atomic beams, W = B/a^, where a is the periodicity of the 
diffraction grating, and B ~ 0.5 at room temperature IHEl- The rule of thumb 02] is that 1^/12 < 0.1, for 
sharp, elastic diffraction patterns to be seen. The requirement for observing interference is more stringent, 
however; the probability of emitting even a single phonon during the diffraction process must be negligible, 
and W/12 < 0.01 is required |l67j. For nanofabricated diffraction grating and mirrors, 1^/12 ~ lO-'^ at room 
temperature and decreases at lower temperatures. Thus it is highly probable that the zero-phonon process 
is the dominant one, and therefore quantum phase coherence is expected in the interferometer. 

As stringent as this zero-phonon-emission condition on atom interferometers may be, atom interferome- 
ters of the same type as MIGO have been successfully constructed before. For example, the interferometer 
used by Pritchard, which is very similar to the horizontal MIGO configurations, uses 1027 K sodium atom 
beams using 100 nm nanofabricated diffraction gratings. Both the temperature of the atoms and the size 
of the diffraction gratings are comparable to MIGO's. The fact that fringes have already been observed in 
Pritchard-type interferometers shows that decoherence due to the atom-beam-splitter interaction does not 
prevent interference from happening. 

4.5 Measuring New General Relativistic Effects 

The underlying advantages of MIGO over LIGO or LISA go beyond being smaller and more sensitive. It is 
also possible to explore with MIGO general-relativistic effects that are not accessible with LIGO. 

4.5.1 Measurements of Stationary Riemann Curvature 

As we have outlined in the Introduction, and have shown in figure ^b), the parallel transport of a four- 
vector around a closed path is a measure of the local Riemann curvature tensor of the spacetime. Referring 
to figureEl E^nd extending the n = ±1 paths that an atom can take into spacetime, it becomes clear that in 
MIGO the atom does precisely what is indicated in figure ^b): It is parallel transported around a closed 
path in spacetime, and as a result its phase is shifted by the local curvature of the spacetime j42| through an 
anholonomy. Although the delta-function impulsive forces in the mirrors impart an additional acceleration- 
and time-dependent phase shift to the atoms, the Riemann-curvature-tensor contribution still remains, and 
will be dominant for low-frequency GWs, and for stationary Riemann curvature tensors in general. It is thus 
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not surprising that MIGO is sensitive to the total local Riemann curvature tensor of spacetime, independent 
of its source, and not just the fluctuations caused by the passage of a GW, as LIGO is. Indeed, while in our 
analysis in Appendix B we have for clarity taken the acceleration due to gravity to be a constant — thus 
neglecting the contribution of the curvature from stationary sources such as the Earth — if we repeat the 
analysis focusing now for stationary gravitational fields we find that the horizontal MIGO configuration will 
measure a phase shift 021 

\^K%o\ = ^^^±^11^^*0,, (13) 

where Rf^^oy the total local Riemann curvature tensor from stationary sources (once again the Sagnac 
effect is neglected). The Earth-based horizontal MIGO can thus measure RoxOy > 4 x 10~^'* s~^ in a one 
second integration time. With a local curvature of 1.23 x 10~^ s~^, it will certainly be possible to measure 
the Earth's Riemann curvature tensor if the contribution from the Sagnac effect can be subtracted from 
the signal independent of the Earth's curvature's. It will also be able to measure the Riemann curvature 
tensor of the Sun (|i?oiOj | = 3.2 x IQ-^'' s'^), and the Moon {\RoiOj \ = 7.0 x 10"" s'^) with MIGO; both 
curvatures vary predictably with time, and their signal can be discriminated from the static phase shifts 
cause by the Sagnac effect and the Earth. For comparison, a GW with an strain amplitude of 10"^^, and a 
frequency of 10"' Hz from a supernova source shown in figuredhas a rms |i?ojOj| ~ 2.5 x lO^^** s^^, similar 
to the nearly static Riemann curvature of the Sun and Moon. If it is possible to subtract out both the 
Sagnac and Earth's curvature contributions, at a |i?oiOj| ^ 1 x 10"^^ s~^, the Earth-based horizontal MIGO 
will be able to measure the Lense-Thirring field of the Earth as well. 

On a classical level, this Riemann-curvature-driven effect on the atom is surprising. In general relativity, 
an atom is modelled as a spin-zero, point-like test particle, and thus should not be directly affected by the 
curvature. However, on the quantum level the atom is described by a delocalized wavefunction ip, which is a 
section of a U{1) line bundle over the local spacetime manifold. Curvature-driven effects of the U{1) fibers 
of the line-bundle above a flat base manifold is well known in quantum mechanics, and is the fundamental 
cause of Berry's phase jl5l I17j . It should not be surprising that the Riemann curvature of a curved base 
manifold will cause phase shifts as well, in a gravitational version of the Aharonov-Bohm effect |12j . 

4.5.2 Spatial Variations of GWs 

The length of LIGO's arms are such that Ll/go ^ Amm/27r, where Xmin is the shortest wavelength of GWs 
in LIGO's frequency range. Variations in the GW in space can then be neglected, and only variations in 
time matter, leading to Eq. Q. However, if the interferometer is larger than wavelength of the GW, general 
relativity predicts that additional terms that depend on the velocity of the particle as well as the spatial 
variation of the GW — the magnetic part of the Weyl tensor of the GW — will now appear in Eq. . It is in 
principle possible to measure these new effects with MIGO; at 17 km, the length of the Space-based MIGO 
is comparable to the wavelength of a GW with a frequency of roughly 3 kHz. As can be seen in figured 
there are a number of GW sources above this frequency that can potentially be used to probe these effects. 

5 Conclusions 

In the first part of this paper, we have explored the implications, broadly speaking, of non-Euclidean 
geometry for quantum physics. The problem of the interaction of gravitational radiation with quantum 
matter, including nonrelativistic quantum many-body-systems, was explored. In the second part of this 
paper, we have explored in detail one particular aspect of this problem, namely, the detection of astrophysical 
sources of the GWs by means of matter- wave interferometry (MIGO). Our theoretical analysis indicates 
that in terms of sensitivity, size, capability, and flexibility, the quantum methods embodied in MIGO have 
overwhelming advantages over the classical methods embodied in LIGO in studying general relativistic 
effects. What before took many kilometers can now be done in a couple hundred meters, and what took 
millions of kilometers can now be done with only thousands of meters. The impact of MIGO will go beyond 
simply being another astronomical instrument, however. With the advantage in MIGO that slowly moving 
atoms are predicted to have over ultrarelativistic ones, it will also spur theoretical and experimental studies 
that probe the intersection of nonrelativistic quantum mechanics with general relativity |17| . As an example, 
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we note that limitations to MIGO's sensitivity are, at their most fundamental level, imposed by quantum 
mechanics, and thus cannot be exceeded. Since MIGO in its essence measures local curvature, like Wigner 
|68[ 169) . we ask: What are the fundamental quantum limitations to the measurement of the Riemann 
curvature tensor arising from the uncertainty principle? This leads naturally to the age-old question, 
approached now from a different angle: What does it imply for Einstein's theory of general relativity, and 
for geometry as a whole, if curvature cannot be measured beyond these fundamental quantum limits? This 
question, which, in order to answer, was once thought to require energy scales on the order of the Planck 
energy — the energy where general relativity and the other forces of nature are expected unify — can now be 
addressed using matter- wave interferometry, once the fundamental quantum limits to MIGO's sensitivity 
have been established. 
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A Brief Review of Gravitational Waves 

In this appendix, we present a brief review of GWs in the linearized gravity limit; the reader is referred to 
^ or UOj for a complete review. 

Given a spacetime manifold M, the metric 5^,^ gives the measure of length, ds — 5/^1/ 

(Ixf^dx", on M. 

The precise form, of g^^, depends on the coordinate system chosen, and if a;^ is another choice of coordinates, 
then the coordinates a;^ expressed in the coordinates is simply Xfj.{x). Since the length of a vector in Ai 
cannot depend on the choice of coordinate systems, g^i,{x)dx^dx^ = g ^v{x)dx^ dx^ , and by using the chain 
rule, we see that in a coordinate transformation the metric changes by 

dx^ dx^ 

g^^ii)=g^,i:,) — — . (14) 

No physically measurable quantities can depend on the choice of coordinates. In Newtonian gravity, and 
in special relativity, we require physics to be invariant under global coordinate transformations. In general 
relativity, we require the theory to be invariant under general coordinate transformations. Consequently, 
the choice of coordinate systems is often called a 'gauge choice', in analogy to EM, and general coordinate- 
transformation-invariance is often referred to as 'gauge invariance' 

We are interested in the propagation of GWs in the spacetime surrounding the Earth. Because of the 
relatively small mass of the Earth, the spacetime can be approximated as being flat. If 77^1, = diag{—l, 1, 1, 1) 
is the Minkowski metric, we can approximate g^i, « r]^i, + /i^j^, where the components of h^i, are 'small' 
compared with unity, and h^i, represents the 'strains' in the flat spacetime caused by the passage of a GW. 
It is treated as just another tensor on A4, and its indices are raised and lowered using 77^1,: h"^ — rj'^^hafj,- 
In linearized gravity, we keep terms only linear in /i^^; since g^"gav = S'^: g^^ ~ V^^'^ — h^"^ . 

Even though we are working in the linearized gravity limit, there is still a remnant of the coordinate 
transformation invariance Eq. H14() left. If x^ ~ x^ + where is a 'small', arbitrary vector, then from 
Eq. (|14|l . hfj^i, transforms to h^j^^ = h^^ + d^£,u + du£,^. The similarity between the transformation of /i^^ under 
a coordinate transformation, and the transformation of the vector potential —^ + d^(f> under a U{1) 
gauge transformation is readily apparent. Being a symmetric, 2x2 tensor, h^^, has ten components. We have 
the freedom, as with the vector potential in EM, to make specific choices for and its derivatives — eight 
variables in all — that will determine our coordinate system. Thus, /i^j^ contains only 10 — 8 = 2 physical 
degrees of freedom, or polarization states. The usual coordinate choice is the TT gauge: rj^^'^h^^i, — 0, 
= 0, and, for GWs in Minkowski space, the condition /iq^ = is automatically follows. 
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In the linearized gravity limit, the Levi-Civita connection is in general, 

r^. = \ {d^K + s./i^ - d'^h^.) , (15) 

where 9^ = djdx^ , while the Riemann curvature tensor and the Ricci tensor are, respectively 

were R^^ = Consequently, R = f]^"' R^,u = -D^J^ + d^d^h'"' . For GWs in the TT gauge, the above 

equations reduce even further, and we see that i? = 0, i?o/j = Oi and i?oiOj = ~^y /2- 
The evolution equation for GWs comes from Einstein's equation, 

Rfj,iy - -;^9tiyR = —^T^^, (17) 

where G is Newton's constant, and T^j^ is the energy-momentum tensor of matter in the spacetime. Away 
from all sources, T^i, = 0, and since i? = 0, we find that Rp.i, = 0, or from Eq. H16|l . d'^dahij = 0, with the 
additional transversality condition d^hij = 0. Thus, like an EM wave in the radiation gauge, hij in the TT 
gauge is a transverse wave that satisfies the wave equation. It also can be expanded in plane waves, 

f'^^- E E 4fik)hsik)e^^'-^~-'\ (18) 

s=+,x k^kyk^ 

were (\f (k) are the unit polarization vectors, and the transversality condition becomes k^tf^ {k) ~ 0. Thus, 
in the plane perpendicular to k, elj can be represented by 2 x 2 matrices, where 

Sketches of the force- field lines produced by these polarizations are drawn in figure |3| 



B Calculation of MIGO Phase Shifts 

We consider first the horizontal MIGO configuration figure |5l Assuming that the atoms in figure |S1 are fast 
moving, in the absence of a GW the atoms travel along the straight-line paths as drawn. At time t^, the 
rth atom is diffracted by the beam splitter, and for the n = -1-1 order, the path that the atom travels is: 
VQy = w±, yo(i) = v±(t — tr), and vqx = v\\, xo{t) = v\\{t — tr), for < t < tr + T/2, where T/2 is the time for 
the atom to travel from the beam splitter to the mirror. For tr + T/2 < t < tr + T the path is: voy = —v±, 
yo{t) = L± — v±{t — tr), while vqx = W|| and xo{t) = — tr) still. Clearly, L± — v±T and i|| = v\fT. The 
velocities and ?;|| are the initial velocities of the atom as it leaves the diffraction-grating beam splitter, 
and just like for light, v±_ ~ 27rh/ma where a is the periodicity of the grating. Then V|| = (v^ — v^)^/^ where 
Vs is the velocity of the beam incident the beam splitter. Gravity can be neglected as long as W||T >> gT'^/2 
where g is the acceleration due to Earth's gravity. Since our focus is on GWs, we approximate 17 as a 
constant, and neglect local curvature effects from stationary sources such as the Earth. A similar set of 
equations hold for the n = —1 order with v± —v±. 

If a GW with hij{t) is present, the paths of the atom will be slightly perturbed from straight lines, 
and we take x ^ xq + xi and y = j/o + where xi and yi are deviations from the unperturbed paths Xq 
and j/o- These perturbations satisfy the geodesic deviation equations of motion in Eq. but with {x,y) 
approximated as {xo,yo) on the right hand side. 

To take into account the action of the GW on the mirrors, we model the mirror and its connection to 
the frame of the interferometer as a spring with a resonance frequency /o — which depends on the size of 
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the interferometer as well as its material properties — and a quality factor Q. We denote by {X{t),Y{t)) the 
position of the mirror along the n = +1 path, and since in the absence of a GW the mirror does not move, 
we take X{t) = i||/2 + Xi{t), and Y{t) = + Yi{t)^ where Xi(t) and Yi{t) are perturbations of the 

mirror's position due to the GW. They are also solutions of Eq. (j^l, but because the mirrors are connected 
to springs, 



which are driven simple harmonic oscillators. 

The action corresponding to the geodesic deviation equation is [71] 

j dti^-mv^ - -mhijv'x^ j , (21) 
Writing v = vq + vi, wc find that for the n = +1 path the change in the action due to the GW is 

« L|| ivi^iU + T) - + T) + + T) - [h^yiU + T) - h^yiU + T/2)] j + 



L 



2- [?^iy(ir + T-/2) - viyiU + T+/2)] + ^v^L^hyyiU + T/2)- 

tr + T ( j-tr + Tjl ft-r+T 



{v\Kx{t) + V\hyy{t)^ dt - VJ_V\\ I / h^y{t) dt' Ky{t) dt S , (22) 



tr + T/2 



after successive integration by parts of Eq. H21(l . and using Eq. Q. The velocities viy{tr + /2) and 
viy{tr + T'^ /2) are the perturbed velocities of the atom right before and right after the mirrors, respectively. 
Since the passage of the GW will shift the velocity of the mirrors as well as the atoms, a jump condition of 
the y-velocity at the mirror requires that 

viyiU + T+/2) = -viyiU + T-/2) + ViyiU + T/2). (23) 

Taking the Fourier transform of /i^ , and neglecting the transients in the solution of Eq, H2U|) , we find for a 
horizontal MIGO that 



A<?Igo(/) - --^i±i||*Ax(/)e-*"^^F^(/T), (24) 



for a GW with amplitude hx{f) and frequency /. Only the x polarization causes a net phase shift in the 
interferometer; the + polarization does not, as can be seen due to the inherent bilateral symmetry of the 
interferometer. The horizontal resonance function is 



TTfT 



sine 



TTfT 



2 



1 

^ (25) 



2/2-/o' + i/o//Q' 

where sinc(x) = sin.T/a;. The analysis of the vertical configuration MIGO follows in the same way. The 
overall form of Eq. 124|) still holds, with only the replacement Fh{fT /A) by 2Fy{fT), where now 



E,{fT) = 1- 



smc 



tt/T 



2 



[1 - sinc(^/T)] - i , (26) 



nfT/2' ' ' " 2/2-/2 + z/o//Q' 

for this configuration. 

If we consider the low frequency / ^ limit of Eqs. (|25ll and H26I) . we find that Fii{fT) « —infT, while 
Fy z7r/T/3, irrespective of how the mirror is mounted onto the frame of the interferometer. This is not 
the case in the high frequency limit. We find that in the infinitely-rigid interferometer limit where /o ^ oo, 
and then / — > oo, both resonance factors approach one. If, however, the mirrors are freely suspended like 
those of LIGO, then /o — > first. Then when the / — s- oo limit is taken, both resonance functions approach 
one- half. In the MIGO graphs of figure [TJ we took as the speed of sound for the material used to construct 
the interferometer frame as 6000 m/s, and that the motion of the mirrors were critically damped with Q = 1. 
This resulted in the small dips at the high-frequency end of the MIGO graphs. 
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C Effects of Thermal Noise on MIGO 



We consider the mirrors to be rigidly attached to the frame of the interferometer. (This analysis can be 
repeated with the mirrors attached to springs as described in Appendix B). The finite temperature of 
the mirror induces fluctuations on the surface of the mirror, i. e., transverse surface sound waves, which 
will induce random Doppler-like shifts in the velocity of the atoms after it hits the mirror. Such velocity 
shifts are not important for LIGO because the speed of light is much larger compared to the velocity of 
these fluctuations. For slowly moving atoms, on the other hand, these fluctuations are important, and can 
potentially degrade the sensitivity of MIGO. 

We model the surface of the mirror as a continuous sheet with transverse surface fluctuations u{x,t) = 
Mosin(fcs • X — (jJst + 4>s), where x = {x,y) is a vector at the surface of the mirror, is the wavenumber of 
the sound mode, Us is its frequency, and 0s is its random phase. The dispersion relation for the sound wave 
is ujs = Vsoundks whcrc Vsound IS the transverse speed of sound for the mirror. The surface of the mirror in 
the presence of the sound mode is defined hy z = u{x, t). 

Consider an atom with wavefunction ipindx, z,t) = AQexp{i{k±z + uj±t)} normally incident on the 
mirror. The reflected wavefunction iprefix, z, t) is determined by the Schrodinger equation, and the boundary 
conditions that the total wavefunction vanishes at z = u(x,t). The total current normal to the surface of 
the mirror vanishes there as well. Because this surface is changing with time, this is a moving boundary 
problem that cannot be solved exactly. We note, however, that for ujs » uj±, the surface of the mirror 
oscillates many times during the time when the atom reflects off of it; these oscillations cancel each other, 
and do not adversely affect the motion of the atom. On the other hand, when uig << i^±, the surface is 
frozen in time, and the problem reduces to an atom reflecting off a rough surface. It is straightforward to 
see that for this case any variations in the reflected atom's velocity normal to the mirror is due to changes 
in the direction of the normal vector to the surface of the mirror. This causes a velocity shift of the atom 
that is proportional to (kgUo)'^, and as we shall see, is not the dominant effect. 

Consider now the case when ujs ~ Expanding V'mc at the surface of the mirror, we get 

'il>ref[x,z^u{x,t),t] = - Aq exp{i {k±u{x , t) + UJ±t)} , 

q—oo 

= -Aq ^ Jqik±uo) exp |i (uj± - q k± ■ x ~ Ugt + (f)s ^ | . (27) 

q—~oo 

The mirror surface acts as a FM modulator, and produces sidebands with an initial amplitude given by the 
gth order Bessel function Jg(fcj^uo). The amplitude of the mode uq — Vth/i-^s, where vth = {kBTmirr/MY^^ 
with Tmirr the temperature of the mirror with mass M, while k± = 27r/a, where a is the periodicity of the 
diffraction grating used in the beam splitter. The product k±uo = k±^vth/^s ~ k±vth/^± = vth/v± << 1. 
Consequently, only the q — 1 sideband is important (since u!± ~ ujg, the q — —l sideband is not a propagating 
mode). 

Thus, since Jo{k±uo) « 1, the incident wavefunction is reflected backward with almost the same am- 
plitude it came in with, but now an additional q = I component with amplitude Ji{k±UQ) mixed in. But 
because the phase of the sound mode is arbitrary, the wavenumber of this reflected wave is approximately 
— fc_L + ki{vth/v±)'^ /2, where fci = {imuij^/hy^'^ is the wavenumber of the q = 1 reflected wavefunction. 
The linear term in the expansion of Ji(fcj^uo) vanishes since roughly half the time the sound mode is in 
phase with the incident atom, and half the time it is out of phase. Thus, the fractional change in the atom's 
velocity is due to thermal motion of the mirror is Av/v± = {vth/v±)'^, and is independent of the material 
properties of the mirror, as expected. This shift in the atom's velocity introduces a thermal-noise-limited 
sensitivity 

^(Z)"^ = 2^A|^^(/T) 1/3/2' ^2^) 

which dies out faster than the shot-noise-limited sensitivity /i(/)s/iot ^ ^/f for high frequencies, but it 
is dominant for low frequencies. The thermal- and shot-noise-limited sensitivities are equal at ftum = 
N {mL ±Av / h)'^ . For a one-kilogram-mass mirror at a temperature of 1 K with a 6000 m/s transverse speed 
of sound, this is roughly 10"^ — 10~® Hz for the ^^^Cs atom MIGOs shown in figure[71 and causes the curves 
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of the Earth-based vertical MIGO and the Space-based MIGO to angle slightly upwards at low frequencies. 
For the Earth-based horizontal MIGO this frequency is roughly 10~^^, and thermal effects are negligible. 
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